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ABSTRACT 



A detailed analysis of more than 800 electronic high-resolution spectra of gamma Cas, which were obtained during a time interval of 
over 6000 days (16.84 yrs) at several observatories, documents the smooth variations in the density and/or extend of its circumstellar 
envelope. We found a clear anticorrelation between the peak intensity and FWHM of the Ho- emission, which seems to agree with 
recent models of such emission lines. The main result of this study is a confirmation of the binary nature of the object, determination 
of a rehable lineal- ephemeris TminRv = HJD (2452081.9 + 0.6) + (203'!52±ff!08)x£, and a rather definitive set of orbital elements. We 
clearly demonstrated that the orbit is circular within the limits of accuracy of our measurements and has a semi-amplitude of radial- 
velocity curve of 4.30±0.09 km s"'. No trace of the low-mass secondary was found. The time distribution of our spectra does not 
allow a reliable investigation of rapid spectral variations, which are undoubtedly present in the spectra. We postpone this investigation 
for a future study, based on series of dedicated whole-night spectral observations. 

Key words, stars: early-type - stars: binaries - stars: Be - stars: individual: gamma Cas 



1 . Introduction 

The well-known Be star of spectral type BOIVe y Cassiopeiae 
(27 Cas, HR 264, HD 5394, HIP 4427, MWC 9, ADS 782A), is 
one of the first two Be stars ever discovered (see ISecchi|[T866l) 
and a member of a visual multiple system. It exhibits spectral 
and brightness variations on several timescales. It underwent 
two major shell phases in 1935-36 and 1939^0. Afterwards, 
it appeared briefly as a normal B star. Emission strength of the 
Balmer lines and the brightness of the star in the visual region 
had been rising slowly during the rest of the 20* century. The 
ob servational history of the star has been summarized in detail 
bv lHarmane3(l2002l) . 

In 1976, y Cas was identified as an X-ray source. This dis- 
covery started a long debate over whether the source of X-rays 
is the star itself or whether y Cas is an X-ray binary with a mass 
accreting compact binary companion. I n an effort to prove the 
duplicity of y Cas, ICowIev et alJ (Il976 f) measured radial veloci- 
ties (RVs hereafter) on a rich collection of photographic spectra 
obtained in the years 1941-1967. They were unable to find any 
RV changes exceeding 20 km s ' or to detect any coherent peri- 



ods between 2'?5 and 400a?0. lJaradl(ll987l) measured RVs on 8 1 
medium-dispersion (30 A mm"') photographic spectra using the 
cross-correlation technique. They concluded that the RVs vary 
with a short period of 1^16885, a semi-amplitude of 27.7 km s"', 
and a well-define d phase curve. Com bining their RVs with 
those measured by ICowIev et al.l (Il976 l). they found a period of 
O?705163 with a semi-amplitude of only 8.6 km s They pre- 
ferred the shorter period, which they interpreted as either arota- 
tional or pulsation period of the star. Robinson & SmithI (l2000h 
published a detailed study of the X-ray flux of y Cas. They found 
that the X-ray flux varied with a period P = 1^12277, which they 
tentatively interpreted as the rotational period of y Cas. They 
used this finding as one of the arguments against the binary 
scenario for the X-ray flux. More recently, ISmith et aTl (l2006h 
have reported a coherent periodicity of l'?21581+0'?00004 from 
the 1998-2006 optical photometry, prewhitened for variations on 
longer time scales. The latter authors pointed out that the initial 
period estimate of 1^12277 was probably an alias of the correct 
period near 1^21581. 
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iHarmanec et al.l (1200 0*) measured RVs of the steep Ha emis- 
sion wings in a series of 295 Ondfejov Reticon spectra spanning 
nearly 2500 days from 1993 to 2000. After removing the long- 
term RV changes, they discovered periodic RV variations with a 
period P = 203'?59+O?29, semi-ampHtude Ki = 4.68 km s"', and 
eccentricity e = 0.26, which they interpreted as the binary motion 
around a common centre of gravity with a low-mass companion. 
They demonstrate that the published RVs from the photographic 
spectra can also be reconciled with the 203^59 period and discuss 
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th e possible properties of the sys tem. Their result was confirmed 
bv lMiroshnichenko et alJ (l2002h . who also measured the RVs of 
the Hot emission wings in a series of 130 electronic echelle spec- 
tra, secured with the 1 m reflector of the Ritter Observatory be- 
tween 1993 and 2002. These two studies differ in the technique 
of RV measurements. While Harmanec et al. (2000) measured 
the RVs manually, sliding the direct and reversed continuum- 
normalized line profiles within a range of intensities on the com- 
puter screen until the best match was obtained, Miroshnichenko 
et al. (2002) also matched the original and reversed profiles, 
but used an automatic procedure. They arrived at a period of 
205'?50+O?38, semi-ampHtude of 3.80+0.12 km s"', and a cir- 
cular orbit. They discuss several possi ble reasons why their re - 
sults differ significant l y frorn those of iHarmanec et aL I (l2000l) . 
iMiroshnichenko et alj (|2002 |) also document the cyclic long- 
term spectral variations of y Cas over a time interval from 1973 
to 2002. 

To shed more light on the differences between these two RVs 
studies, to obtain truly reliable orbital elements of y Cas, and to 
exclude possible 1 d aliases of the 204 d period, we combined 
our efforts and analysed the two sets of spectra, complemented 
by more recent observations from Ondfejov and additional spec- 
tra from the Dominion Astrophysical (DAO), Haute Provence 
(OHP) and Castanet-Tolosan Observatories. The RVs in all these 
spectra were measured by both measuring techniques - alterna- 
tively used by Harmanec et al. (2000) and Miroshnichenko et al. 
( I2OO2I) - and analysed in several different ways. Here we present 
the results of our investigation. We also studied the long-term 
and rapid spectral variations of 7 Cas in our data. 

2. Spectral variations 

We have collected and analysed series of electronic spectra from 
five observatories. They cover the time interval from 1993 to 
2010. The journal of the observations is in Table [1] where the 
wavelength range, time interval covered, the number of spectra, 
and the spectral lines are given. For more details on the individ- 
ual datasets, readers are referred to AppendixlAl 



Table 1. Journal of spectral observations. 



Origin of 




AT 


Lines 


N 


spectra 


(A) 


(RJD) 






Ond 


6300-6730 


49279-55398 


Ha, He, Si 


439 


Rit 


6528-6595 


49272-52671 


Ha 


204 


OHP 


4000-6800 


50414-52871 


Ha, He, Si 


34 


DAO 


6155-6755 


52439-54911 


Ha, He, Si 


136 


OHP/Cst 


various 


53997-55422 


Ha 


13 



Notes. Ond = 2 m reflector of the Astronomical Institute AS 
CR Ondrejov, Rit = 1 m reflector of the Ritter Observatory of 
the University of Toledo, DAO = 1.22 m reflector of the Dominion 
Astrophysical Observatory, OHP = 1.52 m reflector of the Haute 
Provence Observatory, Cast = Castanet-Tolosan; AA = the wavelength 
region covered, AT = the time interval spanned by each dataset, where 
times are given in the reduced Julian dates RJD = HJD-2400000, Lines: 
Ha = Ha, He = He I 6678 A, Si = Si II 6347 A, and Si II 637 1 A. 

We focused our study on the lines in the Ha region, which 
are available for all spectra, although several echelle spectra 
cover almost the whole visible region of the electromagnetic 
spectrum. In particular, we studied the following spectral lines: 
Ha, Hex 6678 A, Sill 6347 A, and Sill 6371 A. No dramatic 



changes were found in these line profiles. The He I 6678 A and 
Si II lines exhibit double-peaked emissions with the well-known 
V/R variations (changes in the relative strength of the shorter 
wavelength, "violet", to the longer wavelength, "red", peak) on 
the timescale of several years. Over the whole time interval cov- 
ered by our spectra, the Ha line was observed as a strong, basi- 
cally single-peaked emission, having a peak intensity between 
3.5 and 5.0 in the units of the continuum level. Its V/R varia- 
tions manifest themselves as a relative shift in the emission peak 
with respect to the centre of the emission profile. Several shal- 
low absorptions can be noted in the Ho- line in some of the stud- 
ied spectra, but most of them are the telluric water vapour lines. 
Ocassionally, some weak shallow absorptions of probably stellar 
origin were seen, but they disappeared in less than several tens 
of days, and we found no regularity in their appearance and dis- 
appearance. The He I 6678 A line consists of a double-peaked 
emission filling a large part of the rotationally broadened pho- 
tospheric (or pseudophotospheric) absorption. The whole line is 
very weak and can only be measured reliably on the spectra with 
high S IN. The emission peaks rise only a few percent above the 
continuum level. Nevertheless, the time variations are seen most 
prominently in this line. The red peak of the He I 6678 A emis- 
sion disappeared almost completely at certain times. It is hard 
to say whether these variations represent only real long-term 
changes or whether they are partly caused by line blending. The 
Sill 6347 A and Sill 6371 A double emission lines are even 
weaker than the He I 6678 A Une, and their peak intensity never 
exceeds 5% of the continuum level. The more recent evolu- 
tion of the Ha line profile is shown in Fig. [T] All line profiles 
shown were ob tained after RJD = 522 25 and were not included 
in the study bv IHarmanec et all (l2000i) . A similar sequences of 
the He I 6678 A, Sill 6347 A, and Sill 6371 A Une profiles 
are shown in Fig. |2] All displayed spectra are from Ondfejov, to 
compare the data with the same resolution. There are, however, 
the huge differences in the flux scale between Ha, He I 6678 A, 
Si II 6347 A, and Si II 6371 A fines. 

8 I < 1 1 < , 1 
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- 

_J I , 1 1 , , 

6532 6542 6552 6562 6572 6582 6592 
wavelength[A] 

Fig. 1. Recent evolution of the Ha line profile. Mid-exposure 
times of the displayed spectra are in RJD = HJD-2400000. 

We fitted the Ha line profiles with a Gaussian profile to 
obtain their peak intensity (/p hereafter) and full width at half 
maximum (FWHM hereafter). This procedure naturally returns 
a value of /p, which is slightly less than the very maximum of 
the emission profile, which is affected by both the VjR varia- 
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6667 6677 6687 
wavelength[A] 



6331 6341 6351 6361 6371 6381 
wavelength[A] 



Fig. 2. Recent evolution of the He I 6678 A (left). Si II 6347 A, 
and Sill 6371 A (right) line profiles. Mid-exposure times of the 
displayed spectra are in RJD = HJD-2400000. The vertical axis 
is in the units of the continuum flux. 



tions and blending with the neighbouring telluric lines. We do 
believe, however, that the fitted Gaussian provides an objective 
measure of the gradual changes in the emission-line strength. We 
omitted the saturated Ha line profiles, of course. One example 
of a Gaussian fit is in Fig. |3]to show where the FWHM and /p 
were measured. 



continuum level 




6560 6570 
waveIength[A] 



Fig. 3. An example of the Gaussian fit to an Ha line profile, 
which also shows the derived quantities /p and FWHM. Dashed 
line: the observed Ha profile; solid line: the Gaussian fit to it; 
dotted line: the continuum level. 



Figure |4] shows the time variations of the FWHM and /p. 
An interesting finding is that the secular variations of these two 
quantities are anticorrelated with each other The apparently in- 
creased scatter of both dependencies between RJDs w 50000 and 
52000 is caused solely by the lower resolution in intensity of the 
Ondfejov Reticon spectra taken prior RJD = 52000.0 



' Although the original Reticon detector and the currently used CCD 
detector were attached to the same camera of the coude spectrograph 
and have the same spectral resolution, the control electronics of the 
Reticon detector allowed distinguishing only 4000 steps in intensity, 
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Fig. 4. The secular time variations of /p (upper panel) and 
FWHM (bottom panel) of the Ha line. The time on abscissa is in 
RJD = HJD-2400000. The open symbols denote measurements 
from the Ondfejov Reticon detector, capable of distinguishing 
only 4000 intensity steps, which is why these measurements de- 
viate systematically from the rest. See the text for details. 
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Fig. 5. Comparison of the automatically measured Ha emission- 
wing RVs (ordinate = y) with those measured manually (abscissa 
= x). The dashed fine is a fitted linear function y - 0.968 x. 



3. RV changes 

3. 1. Long-term and periodic RVclianges and tlie new 
orbital solutions 



Simila r to iHarmanec et al.l (l2000t) and iMiroshnichenko et al.l 
(2002'), we measured the steep wings of the Ha emission in all 
unsaturated profiles. The reasons the emission wings should pro- 
vide a good estimate of the true orbital motion of the Be pri- 
mary around the common centre of gravity with the secondary 
were recently summarized in detail by Ruzdiak et al. (2009). To 
them, we can add that Poeckert & Marlbor ough! (i 1 978 ) modelled 
the Ha emission of y Cas, and their model showed that the Ha 
emission wings originate in regions that are much closer to the 
star than the radiation that is forming the upper part of the line. 



while the CCD detector recognizes 60000 intensity steps. This leads to 
a systematic difference for strong emission lines. 
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Fig. 6. Plots of RVs vs. RJD = HJD-2400000: Top panel: 
Manually measured Ha emission-wing RVs. Bottom panel: 
Automatically measured Ha emission-wing RVs. The solid lines 
in all panels represent the long-term RV change as derived with 
the program HEC13. The squares in both panels denote the lo- 
cal 'systemic' velocities calculated with the program SPEL for 
individual data subsets. The dashed line represents Hermite- 
polynomial fit computed with program HEC36. See the text for 
details. 

Moreover, we also attempted to measure the RVs of 
other available specti-al lines (He I 6678 A, Si II 6347 A and 
Sill 6371 A) to see if they undergo similar time changes. We 
primarily measured also the outer emission wings of these lines, 
but for the He I 6678 A line it was possible to obtain a relatively 
accurate RV of the central absorption core. The RVs measured 
on the emission wings of the Si II 6347 A and Si II 6371 A lines 
were averaged. 

Two methods of RV measurement were used as follows. 

1. Manual measurements were carried out in the program 
SPEFO, written by Dr. J. Horn and more recently improved 
by Dr R Skoda and Mr. J. Kipata (see ' Horn et"anil996t 
[Skodal[T9"96) . The user can slide the direct and reversed line 
profile on the computer screen until a perfect match of the 
selected parts of the profiles is achieved. The advantage of 
this, admittedly a bit tedious procedure, is that the user actu- 
ally sees all measured line profiles and can avoid any flaws 
and blends with the telluric or weak stellar lines. It is invalu- 
able for measurements of weak spectral lines, where any au- 
tomatic method can be easily fooled by the noise. Th e same 
measuring technique has also been used by Harman ec et al.l 
(2000). One of us, JN, also re-measured all Reticon spectra 
used in their study. Plots of the old vs. new measurements are 



available as Fig. IA.4l in Appendix |A] Considering the good 
agreement of both measurements, we used the mean value 
of the original and new RV measurements for each studied 
feature. All lines (whenever available) were measured with 
this technique. 

2. Automatic measurements were obtained with a program 
written by AM, which also shifts the direct and reversed line- 
profile images for a selected range of relative intensities in 
the continuum units to find a minimum difference between 
them. The advantages of this method are that it is fast and im- 
personal. A potential danger is that it can be fooled by flaws 
and blends in some particular cases. Only the Ha emission 
wings were measured with this method. 

We denote the RVs measured by the first method as man- 
ual and those measured by the second method as automatic to 
distinguish them in following sections. In Fig. |5]the automatic 
Ha emission RVs are plotted vs. the manually measured ones 
to see whether there is any systematic difference between the 
two methods. We fitted the data with a linear relation and some- 
what surprisingly the slope was found to be 0.968+0.009; i.e., 
the automatic method finds a slightly narrower total range of RV 
variations than the manual one. All individual RV measurements 
on the steep wings of the Ha emission are published in detail in 
Table 2 for the manual, and in Table 3 for the automatic mea- 
surements Q 

The Ha emission RVs measured manually and automatically 
are plotted vs. time in Fig. |6l Additional time plots for other 
measured features can be found in Appendix|A] Figure |6] shows 
that y Cas exhibits long-term RV variations over several years, 
which seem to correlate with those of the peak intensity of the 
emission. To be able to search for periodic RV changes on a 
shorter timescale, one has first to remove the long-term ones. 
To check how robust the result is or how much it depends on 
the specific way of secular-changes removal, we applied three 
different approaches to this goal. 

The first was to use the program HEC13 written by 
PH , which i s based on the smoothing technique developed 
by IVondrakI (ll969l Il97_7() and which uses some subroutines 
kindly provided by Dr. Vondralfl The level of smoothing is con- 
trolled by a smoothing parameter e (the lower the value of e, the 
higher the smoothing), and the smoothing routine can operate 
either through individual data points or through suitably chosen 
normal points, which are the weighted mean values of the ob- 
served quantity (RV in our case) over the chosen constant time 
intervals. In both cases, the 0-C residua are provided for all 
individual observations. In these particular cases the following 
specifications for smoothing were used: e = 5xlO"'^ and 200 d 
normals for the Ha emission-wing RVs measured by both meth- 
ods, and e = IxlO"' and 200 d normals for the RVs measured 
on the He I 6678 A absorption core. These particular choices 
of £ make the smoothing function follow only the secular RV 
variations. The solid lines in the three panels of Fig. |6] show the 
estimated long-term changes, which were then substracted from 
the original RVs. These prewhitened RVs were than searched for 
periodicity from 3000^0 down to 0^5 with the program HEC27 
(also written by PH), based on the PDM technique developed by 
Stel lingwerj ((1978). The 0-statistics periodograms for the manu- 
ally and automatically measured Ha emission RVs are plotted in 
Figs. |7] and [8] respectively. The upper panels in both plots show 



^ Tables 2 and 3 are published only in the electronic form. 
^ The program HEC13 with brief instructions how to use it is avail- 
able to interested users at http://astro.troja. mff.cum.cz/ftp/hec/HEC13, 
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the range of periods from 3000^ down to 50^, while the pe- 
riods from 2'to down to 0'?5 are shown in the lower panels. The 
periodograms are flat, with 6 close to a value of 1 for all periods 
between 2 and 50 d, which is why we do not show these parts 
of the periodograms in the figures. One can see that the com- 
bination of RVs from several observatories, which are different 
from each other in their local times, safely excluded the one-day 
aliases. 
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Fig. 7. Stellingwerf 's statistics for the manually measured RVs 
of the Ha emission wings plotted vs. frequency /. Upper panel: 
Periods from 3000?0 down to 50!0. Bottom panel: Periods from 
2?0 downto 0^5. The panels have a different scale on the ordinate. 
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Fig. 8. Stellingwerf 's 9 statistics for the automatically measured 
RVs of the Ha emission wings plotted vs. frequency /. Upper 
panel: Periods from 300(y?0 down to 50^0. Bottom panel: Periods 
from 2^0 downto (fiS. The panels have a different scale on the 
ordinate. 



The deepest minimum in both periodograms at a frequency 
/ ^ 0.004910 d"' corresponds to a period ofP^ 203'?0. The two 
shallower peaks at lower frequencies in Figs.|7]and[8]correspond 
to the integer multiples of the 203 d period. The 203 d period 
was also detected in the measured RVs of other spectral lines, 
though with a larger scatter in the RV curves. These additional 
results are presented in the online Appendix. 



Since no obvious signs of the secondary companion are seen 
in the spectrum, we adopted y Cas as a single-line spectro- 
scopic binary. We derived a number of orbital solutions for the 
Ha emission RVs, prewhitened for the long-term changes with 
HEC13. We used the program SPEL (written by Dr. J. Horn and 
never published) for this purpose. The program ha s alrea dy bee n 
used in several previous studies, e.g., Harmanec (1983, 1984]), 
| Koubskv et al.. Cl985.)„Stefl et all ([l990.) . and . Hornet al.. (.19921 
11994. 



Our first goal was to dec i de wh ether the orbital eccentric- 
ity found by Ha rmanec et al.l ( l2000h is r eal or whether the or- 
bit is actually circular as concluded by iMiroshnichenko et al.l 
(l2002h . In Table |4] the eccentric-orbit solutions for the manually 
and a utomatically measured H a emission-wing RVs are com- 
pared. Lucv & SweenevI (Il97lh have pointed out that observa- 
tional uncertainties may cause the estimated eccentricity to be 
biased when the eccentricity is low. The probability that the true 
eccentricity is zero can be calculated, and this is given in the col- 
umn "L-S test". If this probability is greater than 0.05 we accept 
the hypothesis that the true eccentricity is zero at a 5% confi- 
dence level. 

To shed more light on the problem, we split both manually 
and automatically measured RVs into a number of data sub- 
sets, each of them covering a time interval not longer than three 
consecutive orbital periods and containing enough observations 
to define the orbital RV curve. We used the original, not the 
prewhitened RVs. The phase diagrams for a period P = 203^52 
for all selected data subsets are shown in Fig.|9] We derived the 
elliptical-orbit solutions for them, again testing the reality o f the 
non-zero orbital eccentricity after iLucv & SweenevI (Il97lh . To 
always find the solution with the smallest rms error, we started 
the trial solutions for each subset with initial values of oj for 
four possible orientations of the orbit, namely 45°, 135°, 225°, 
and 315°. 

The corresponding orbital elements for the selected subsets, 
together with the Lucy-Sweeney test, are summarized in Table|5] 
for manually measured RVs and in Table |6] for automatically 
measured RVs. The results show very convincingly that the true 
orbit must be circular (or has a very low eccentricity, which is 
beyond the accuracy limit of our data). Although the L-S test 
detected a definite eccentricity for several subsets, the individual 
values of the longitude of periastron are basically accidental, and 
that is probably the strongest argument for the eccentric-orbit so- 
lutions not being trusted. 



Table 4. Eccentric-orbit solutions based on the Ha emission- 
wing RVs measured manually and automatically and 
prewhitened for the long-term changes with the program 
HEC13. 



Solution No.: 


1) 


2) 


Element 


Manual 


Automatic 


F(d) 


203.36±0.10 


203.08±0.11 


T^RVmin (RJD) 


52083+17 


52080±13 


can 


19±30 


147±22 


e 


0.048±0.027 


0.072±0.030 


Ki (km s"') 


3.88±0.10 


3.93+0.12 


y (km s"') 


0.164+0.071 


0.180+0.081 


rms (km s"') 


1.772 


1.948 


L~S test 


0.183 


0.064 


No. of RVs 


757 


700 
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Table 5. The eccentric-orbit solutions for the subsets of manu- 
ally measured Ha emission-wing RVs, shown as phase plots in 
FiglH 





e 


CO 


Ki 


L-S 


A' 






n 


(km s-') 


test 




2 


0.052±0.073 


141+38 


4.30+0.24 


0.734 


25 


3 


0.313±0.082 


165+12 


5.43+0.44 


0.001 


42 


4 


0.416±0.083 


108.9+8.3 


6.46+0.58 


0.000 


42 


5 


0.14±0.14 


160+41 


4.37+0.60 


0.672 


25 


7 


0.10±0.11 


322+47 


6.88+0.88 


0.674 


50 


9 


0.109+0.067 


201+36 


4.10+0.30 


0.299 


42 


10 


0.120+0.040 


263+16 


5.20+0.20 


0.026 


61 


12 


0.086+0.055 


149+28 


4.64+0.26 


0.353 


46 


13 


0.18 ±0.10 


308+29 


3.56+0.36 


0.270 


70 


14 


0.102+0.040 


147+22 


4.31+0.21 


0.038 


59 


16 


0.116+0.067 


105+31 


3.87+0.22 


0.287 


30 



Notes. = a number of RV subset (the same numbering being also 
used in Fig. L-S test = probability that the eccentricity found is 
false, A' = number of RVs in the subset. 



Table 6. The eccentric-orbit solutions for the subsets of auto- 
matically measured Ha emission-wing RVs. 





e 


a> 


Ki 


L~S 


N 






n 


(km s-') 


test 




2 


0.30+0.14 


322+24 


4.57+0.81 


0.174 


24 


4 


0.42+0.11 


285+17 


5.82+0.78 


0.004 


44 


5 


0.34+0.10 


6+20 


4.74+0.62 


0.014 


25 


7 


0.16+0.12 


87+49 


6.0+1.1 


0.311 


52 


9 


0.15+0.19 


92+53 


3.38+0.50 


0.583 


40 


10 


0.098+0.050 


48+23 


4.12+0.19 


0.147 


60 


12 


0.40+0.21 


336.0+9.5 


5.9+1.2 


0.054 


36 


13 


0.21+0.14 


82+44 


3.94+0.53 


0.157 


59 


14 


0.167+0.068 


111+20 


4.57+0.33 


0.026 


46 


16 


0.26+0.13 


227+24 


4.61+0.56 


0.201 


27 



Notes. Ns = a number of a RV subset (the same numbering being also 
used in Fig.|9]l, L-S test = probability that the eccentricity found is false, 
N = number of RVs in the subset. 



Adopting the circular orbit from now on, we attempted 
to remove the long-term RV variations with a different ap- 
proach. Using the manually measured Ha emission-wing RVs, 
we treated each data subset spanning no more than three orbital 
periods (last four subsets) and two orbital periods (remaining 
subsets) as data coming from separate spectrographs, allowing 
SPEL to derive individual 'systemic velocities' for each subset. 
This way, the orbital solution was again free of the long-term 
changes but they were removed in discrete velocity steps. The 
corresponding orbital solution is in Table |7] and individual ve- 
locity levels (y's) are shown in the top panel of Fig. |6]and listed 
in Table lA.2] in Appendix lAl At first sight, this way of removing 
the long-term changes might seem less accurate, but it leads to 
smaller rms error for the solution than the removal via HEC13. 
There are two reasons for that: 

1. The HEC13 program computes normal points from RV sub- 
sets covering constant time intervals, while for y velocities 
the length of subsets was chosen more suitably to our data 
distribution in time. 

2. The HEC13 program computes the normal points as 
weighted means, but the y velocities computed with the 
SPEL program are elements of the Keplerian orbital model. 
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Fig. 9. Phase diagrams for subsets of Ha emission-wing RVs 
measured manually. The ordinates of all plots are RVs in 
km s 'and the individual subsets vertical axis is in RV and the 
different ranges reflect the fact that the original, not prewhitened 
RVs are used. Trial orbital solutions for these subsets are in 
Tabs.|5]and|6l 



For completeness, a solution for an eccentric orbit was also 
derived but the L-S test gave the probability of 0.18, reassuring 
us that the eccentricity is spurious. Also for this method of the 
long-term removal the rms eiTor of the orbital solution is better 
for the manually than for the automatically measured RVs. 

We tested yet another method of removing the long-term 
variations. It can only be used when one has a RV curve uni- 
formly enough covered by observations. The RVs are averaged 
over chosen time intervals, and the Hermite polynomials are fit- 
ted through these averaged (normal) points. We used the pro- 
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Table 7. Orbital elements obtained using RVs measured manu- 
ally (3) and automatically (4) on the Ha line, with the removal of 
the long-term changes using different y velocities for individual 
data subsets. 



Solution No.: 


3) 


4) 


Element 


Manual 


Automatic 


P(d) 


203.65±0.13 


203.47±0.14 


7^RVmin(RJD) 


52081.42+0.81 


52081.99+0.95 


K,{km s-') 


4.084±0.10 


4.14+0.13 


rms(km s"') 


1.657 


1.908 


No. of RVs 


757 


700 



gram HEC23 to compute the normal points and program HEC36 
to fit thenfl We tentatively averaged the RVs over a 300 d and 
a 400 d interval. New orbital solutions were derived using RVs 
prewhitened this way. The rms error of the resulting solution was 
approximately the same as the rms error of the solution for RVs 
prewhitened HEC13. We decided to use this approach in another 
way. We used the systemic velocities derived with SPEL as nor- 
mal points and fitted them with the Hermite polynomials using 
HEC36. The RJDs of RVs in a subset were averaged and the 
mean RJD was used as the epoch of the y velocity. The same 
approach to computing epochs of normal points is also used in 
the program HEC 1 3 . This way we effectively removed one of the 
disadvantages of the previous method since HEC36 connects the 
normal points with a smooth curve, thus removing the discon- 
tinuous shifts introduced with the second method. The Hermite- 
polynomial fit is shown in the first two panels of Fig. |6] The 
0-C residua were again used to derive circular-orbit solutions 
with SPEL. These are presented in Table|8] As expected, the im- 
provement in the resulting rms errors with the second method is 
relatively small. 

Table 8. Orbital elements obtained using the Ho' RVs measured 
manually (5) and automatically (6) and removing the long-term 
RV variations via a Hermite-polynomial fit through the locally 
derived y velocities. 



Solution No.: 


5) 


6) 


Element 


Manual 


Automatic 


P(d) 


203.523+0.076 


203.371+0.089 


rRVmin(RJD) 


52081.89+0.62 


52082.07+0.76 


/TiCkms-') 


4.297+0.090 


4.26+0.11 


7(km s"') 


0.018+0.064 


-0.018+0.075 


rms (km s"') 


1.592 


1.825 


No. of RVs 


757 


700 
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Fig. 10. Top panel: The net orbital RV curve corresponding to the 
circular-orbit solution (5) of Table |8] Bottom panel: The 0-C 
residua from the orbital solution. Both plots have the same ve- 
locity scale, and the epoch of RVmin is used as a reference epoch. 
The short abcissa in the right upper corner of the top panel de- 
notes the rms of 1 observation for the Keplerian fit no. 5. 
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Fig. 11. Top panel: The net orbital RV curve corresponding to the 
circular-orbit solution (6) of Table |8] Bottom panel: The 0-C 
residua from the orbital solution. Both plots have the same ve- 
locity scale, and the epoch of RVmin is used as a reference epoch. 
The short abcissa in the right upper corner of the top panel de- 
notes the rms of 1 observation for the Keplerian fit no. 6. 



The net orbital RV curves and the corresponding 0-C 
residua from the orbital solutions are shown in Fig. [TO] for the 
manual and in Fig. [TT]for the automatically measured Ha RVs. 
The rms errors per one observation of the Keplerian fit no. 5 
based on manually measured RVs (no. 6 in the case of automatic 
measurements) are shown as short abscissae in the upper right 
corners of both figures. The RVs prewhitened for the long-term 
changes, on which these best solutions 5 and 6 are based, are 
also presented in Tables 2 and 3. 

We did several preliminary tests of the possible rapid varia- 
tions of y Cas. Although the results clearly demonstrated their 
presence, the quantitative results were inconclusive so we de- 

* Both programs, written by PH, and the instructions how to use 
them, are available at http://astro.troja.mff.cuni.cz/ftp/hec/HEC36,. 



cided to postpone analysis of rapid spectral changes for a future 
study, based on dedicated series of whole-night spectral obser- 
vations. 

4. Interpretation of results 

4.1. Spectral variations 

We have confirmed the continuing presence of long-term varia- 
tions in the Ha emission profile using the Gaussian fit and RV 
measurements. The former procedure gives more objective re- 
sults than a direct measurement of the peak height and the full 
width of the line at half maximum (FWHM), because the ob- 
served height of the highest peak of a very strong emission line 
partly reflects the long-term V/R changes, while the measured 
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FWHM can be affected by the presence of several telluric wa- 
ter vapour-lines of variable strength. Figure |4] shows that the /p 
and FWHM of the Ho- emission are anti-correlated over the in- 
terval covered by our observations and seem to also show some 
variability on shorter time scales. We carried out period searches 
for both these quantities prewhitened for the long-term changes, 
but no significant periodicities were detected, and trial plots of 
prewhitened data did not show any evidence of variations con- 
nected with the orbital period. 

It is probable that the observed long-term changes reflect 
some variations in the density and/or the extent of the circum- 
stellar disk around the Be prir nary. Using optica l interferometry, 
lOuirrenbach et al.l (Il997h and lTvcner et al.1 (I20b6) were able to 
resolve the envelope around the primary component of 7 Cas and 
estimate that its inclination must be higher than 46° and 55°, re- 
spectively. Adopting a reasonable assumption that the rotational 
axis of the disk is roughly perpendicular to the binary orbit, we 
can conclude that the binary system is also seen under an orbital 
inclination higher than some 45°. We then suggest the following 
possible interpretation of the observed changes. 

It is very probable that the envelope around the Be primary is 
rotationally supported and that its linear rotational velocity de- 
creases with the axial distance from the star Regardless of the 
process causing secular variations of the disk, one can assume 
that the increase in the height of the emission peak reflects the 
presence of more emission power at lower rotational velocities, 
thus implying either an increase in the density of the outer parts 
of the disk and/or an increase in the geometrical extent of the 
disk. This must naturally decrease the observed FWHM of the 
emission as observed. This qualitative scenario seems to be sup- 
ported by the line-profile calculations published by ISilai et al] 
(l2Q10h . 

4.2. Orbital motion 

As mentioned in Sect. [1] one of the principal motivations of this 
study was to resolve the d ifferences in th e orbital solutions ob- 
tained by iHarmanec et al.l (|2000) and by iMiroshnichenko et"an 
(l200l and to arrive at a more definitive set of the orbital el- 
ements. We were able to combine both independent datasets 
and complement them with more recent spectra. In Sect. |3] we 
carried out a number of various tests, analysing separately the 
RVs measured by a manual and an automatic technique. We also 
tested the effect of the different ways of data prewhitening on 
the resulting elements. Special attention was payed to the test of 
whether the orbit is actually circular or has a significant eccen- 
tricity. 

The principal results are the following. 

1. The binary orbit of y Cas is circular, at least within 
the limits of the accuracy of our data, as concluded by 
IMiroshnichenko et aP (l2002h . 

2. The resulting value of the orbital period is now well con- 
strained by the data at hand, and it is robust with respect to 
different ways of analysis. From all experiments, we were 
finding values between 203'?0 and 203'?6 , close to the value 
already found bv iHarmanec et all (l2000l) . An inspection of 
all trial solutions shows that the solutions based on manu- 
ally measured RVs have rms errors that are systematically 
lower for w 15% in comparison to the solutions for the au- 
tomatically measured Ha emission-wing RVs. We therefore 
conclude that solution 5 of Table |8] is the best we can offer 
and suggest the following linear ephemeris for the epoch of 



RV minimum to be used in the future studies of this binary: 

rmm.RV = HJD (2452081.89 + 0.62) 

+ (203'?523 + O?076) X £. (1) 

3. The semiamplitude of the orbital motion is close to 4 km s"' 
for all solutions. The recommended value from solution 5 
is Ki - 4.297+0.090 km s"', implying the mass func- 
tion f(M) - 0.00168 Mq. F or com p arison , Harman ec et al.l 
(2000) and Miroshnic henko et al.1 (l2002h obtained semi- 
ampHtudes of 4.68+0.25 km s'Wd 3.80+0.12 km s"', 
respectively. If we adopt the inclination value i = 45° 
and the primary star mass Mi = 13 Mq suggested by 
[Harmanec et al.l (l2000h , we can estimate the secondary star 
mass M2 = 0.98 Mq. If the system is at a post mass-transfer 
phase, then the secondary might be a hot helium star that 
could be directly detectable in the UV region of the electro- 
magnetic spectrum. 

Note: After our paper was accepted for pubUcation, we had 
the privilege to read a preliminary version of another study of 
-y Cas kindly communicated to us by Dr. Myron A. Smith and 
his coauthors. They analysed in particular a smaller and partly 
independent set of RVs to obtain their own orbital solution. Their 
results are compatible with our final solution within the respec- 
tive error bars. 
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Appendix A: Details of spectroscopic 
observations and tlieir analyses 

A. 1. Observational equipment 

Here we provide more details on the spectra used in this study 
(see Table [T]i and their reduction: 

1. Ondfejov spectra: All 439 spectra were obtained in the 
coude focus of the 2.0 m reflector and have a linear dis- 
persion 17.2 A.mm"' and a 2-pixel resolution R ~ 12600 
(~ 11-12 km s 'per pixel). The first 318 spectra were taken 
with a Reticon 1872RF detector Complete reductions of 
these spectrograms were carried out by Mr Josef Havelka, 
Mr Pavol Habuda and by PH in the program SPEFO. The 
remaining spectra were secured with an SITe-5 800x2000 
CCD detector Their initial reductions (bias substraction, 
flatfielding, extraction of 1-D image and wavelength cali- 
bration) were done by MS with the IRAF program. 

2. DAO spectra: All 136 spectra were obtained in the coude 
focus of the 1 .22 m Dominion Observatory reflector by S Y, 
who carried out initial reductions (bias substraction, flatfield- 
ing, extraction of 1-D image). The wavelength calibration 
of the spectra was carried out by JN in SPEFO. The spec- 
tra were obtained with the 32121H spectrograph with the 
IS32R image sheer. The detectors were UBC-1 4096x200 
CCD for the data prior to May 2005 and SITe-4 4096x2048 
CCD for the data after May 2005. The spectra have a linear 
dispersion of 10 A.mm ' and 2-pixel resolution R ~ 21700 
(~ 7 km s^'per pixel). 

3. OHP spectra: Pubhc ELODIE archiyj^ o f the Haute 
Provence Observatory (iMoultaka et al.l l2004l) contains 35 
echelle spectra obtained at the 1.93 m telescope. They have 
resolution R ~ 42000. Initial reductions (bias substraction, 
flatfielding, extraction of 1-D image, and wavelength cali- 
bration) was carried out at the OHP. We only extracted and 
studied the red parts of the spectra. 

4. Ritter spectra: All 204 spectra were obtained with a fiber- 
fed echelle spectrograph at the 1 m telescope of the Ritter 
Observatory of the University of Toledo. We obtained spec- 
tra in form of ASCII table covering only region close to 
Ha spectra line. The resolution of the spectra is R ~ 26000. 
Initial reductions of spectra (bias substraction, flatfielding, 
extraction of 1-D image and wavelength calibration) were 
carried out at the Ritter Observatory with the IRAF program. 

5. Castanet Tolosan and OHP spectra: We downloaded these 
spectra from the Be Star Spectra databas^ All of them were 
obtained by CB with several different spectrograph^ Only 
spectra with a resolution comparable to spectra obtained at 
the rest of observatories were used in the study. Initial reduc- 
tions (bias substraction, flatfielding, extraction of 1-D im- 
age, and wavelength calibration) were carried out by CB. 

For all individual spectrograms, the zero point of the helio- 
centric wavelength scale was corrected via the RV me asurements 
of sel ected unblended telluric lines in SPEFO (see iHorn et aP 
119961 for details). 



^ URL: http://atlas.obs-hp.fr/elodie/' 
* URL: http://basebe.obspm.fr/basebe/^ 

' For detailed information on the spectrographs used, see the CB 
homepage at http://astrosurf.com/buil/| 



Table A.l. Orbital elements obtained using RVs measured on 
the emission wings of He I 6678 A line, absorption core of the 
He I 6678 A line, and emission wings of the Sill 6347 A and 
Sill 6371 A lines. RJD = HJD-2400000 



spectral 


He I 6678 A 


He I 6678 A 


Sill 


line 


emission 


absorption 


emission 


P(d) 


203.52(fixed) 


203.52(fixed) 


203.52(fixed) 


r„„„(RJD) 


52085.3+3.1 


52081.4+2.6 


52096.4+3.5 


A:i(km s-') 


3.60±0.43 


5.22+0.50 


3.14+0.45 


7(km s"') 


0.16+0.27 


0.18+0.32 


0.29+0.28 


rms(km s"') 


6.179 


7.548 


6.357 


A' 


560 


572 


550 



A.2. Additional RV measurements 

We measured RVs on emission wings and absorption core 
of He I 6678 A and emission wings of Sill 6347 A and 
Sill 6371 A lines. The program SPEFO was used to the task. 
The precision of these RV measurements is quite low, since the 
relative flux in the lines is only several percent greater than in 
surrounding continuum (see Fig. |2]i. One could be easily mis- 
led during measurements, because measured lines are deformed 
with continuum fluctuations, and they blend with telluric lines. 
Despite these complications RVs measured on these lines ex- 
hibit long-term variations very similar to the variations that 
can be seen in Fig. |6] RVs measured on emission wings of 
He I 6678 A line, absorption core of He I 6678 A line, and emis- 
sion wings of Sill 6347 A and Sill 6371 A lines are shown in 
Fig.O 

The long-term variations were removed with the program 
HEC13, using the 200 d normals and e = 5x10 The model 
of long-term variations derived by HEC13 is shown in Fig. lA.ll 
The residua were searched for periodicity using the HEC27 pro- 
gram. A period near 200 d was detected in all cases, although 
with a lower significance than the Ha emission RV (see Figs.|7] 
and[8] The 6 statistics periodograms for trial periods from 3000'!0 
down to 50*^) are shown in Fig. lA.2l where the period P = 203^0 
is denoted. The mininum for this period is not the dominant 
one only in the case of RV measured on the silicon lines. It is 
probably due to their low precision and/or incomplete removal 
of the long-term changes via HEC13. 

The circular-orbit solutions were computed for all 
prewhitened RVs, keeping the orbital period fixed at the 
value P = 203'?52. The corresponding orbital elements are in 
TablelAH 

Data smoothing with different y velocities derived with the 
SPEL program was tested as well. It led to none or a slight (lower 
than 5%) improvement in rms for the Keplerian fit. 

In passing, we wish to mention that the Gaussian fits of the 
Ha line also provided individual RVs of this line. Not surpris- 
ingly, a Keplerian fit of these RVs resulted in a semiamplitude 
lower than was obtained for the directly measured RVs and for 
some 40% greater rms than that for our preferred solution 5. 

A.3. Comparison 

PH and JN measured RVs independently on spectra obtained 
with Reticon detector at Ondfejov observatory. Comparison 
of their results is shown in Fig. IA.4I Dependency was fit- 
ted with linear function y = a.x. The resulting parameters are 
emission wings of Ha line a = 1.024+0.006, emission wings 
of He I 6678 A line 0=1.038+0.024, absorption core of 
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Fig. A.l. A time plot of the RVs measured manually in SPEFO. 
Top panel: The emission wings of the He I 6678 A line, Middle 
panel: The absorption core of the He I 6678 A line. Bottom 
panel: The emission wings of Si II 6347 A and Si II 6371 A lines. 
The HEC13 model of the long-term variations is shown by a 
soUd Une in each panel. 
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Fig.A.2. The 6 statistics periodograms. Top panel: the emis- 
sion wings of He I 6678 A line. Middle panel: the absorption 
core of He I 6678 A line. Bottom panel: the emission wings of 
Sill 6347 A and Sill 6371 A lines. The period ^ = 203^ is 
marked with an arrow. Vertical axis: normalized phase scatter 
6; Horizontal axis: frequency(d '). 

Table A.2. y velocities obtained through the same orbital solu- 
tion as orbital elements in Table Q 



method 


automatic 


manual 


No. of 


7(km s"') 


7(km s~') 


subset 






1 


-16.04±0.29 


-16.54+0.26 


2 


-9.88+0.40 


-9.50+0.34 


3 


-2.54±0.30 


-5.21+0.26 


4 


0.37±0.29 


-1.45+0.26 


5 


0.39±0.39 


-0.70+0.34 


6 


-2.04+0.38 


-1.89+0.36 


7 


-7.36+0.27 


-6.46+0.24 


8 


-9.55+0.36 


-7.89+0.32 


9 


-8.29+0.31 


-9.56+0.26 


10 


-8.81+0.25 


-9.47+0.22 


11 


-9.32+0.43 


-7.89+0.24 


12 


-8.55+0.32 


-9.01+0.25 


13 


-9.25+0.27 


-8.34+0.22 


14 


-7.69+0.29 


-9.12+0.22 


15 


-7.02+0.45 


-9.24+0.39 


16 


-7.43+0.39 


-8.46+0.32 


17 


-7.65+0.53 


-9.51+0.48 



He 1 6678 A hne a = 1 .022+0.008. Differences between the mea- 
surements of both authors are quite high for He I 6678 A line 
emission wings measurements. It is probably because wings of 
the line are affected by the background noise and because the 
red peak of the He I 6678 A line is very low at some point in the 
VI R cycle. 
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Fig. A.3. Phase diagrams of RVs prewithened with HEC13. Top 
panel: the emission wings of the He 1 6678 A line, Middle panel: 
the absorption core of the He I 6678 A line, Bottom panel: the 
emission wings of the Sill 6347 A and Sill 6371 A lines. The 
period P = 203'?52 was used. 
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Fig. A.4. RVs measured by JN plotted vs. RVs measured by PH 

on spectra obtained at Ondfejov Observatory with Reticon detec- 
tor. Top panel: emission wings of Ha line. Middle panel: emis- 
sion wings of He i 6678 A Une, Bottom panel: absorption core 
of He I 6678 A line. 



